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Abstract 

It has been widely claimed that inflation is generically eternal to the future, even in models where 
the inflaton potential monotonically increases away from its minimum. The idea is that quantum 
fluctuations allow the field to jump uphill, thereby continually revitalizing the inflationary process 
in some regions. In this paper we investigate a simple model of this process, pertaining to Xfd 
inflation, in which analytic progress may be made. We calculate several quantities of interest, 
such as the expected number of inflationary efolds, first without and then with various selection 
effects. With no additional weighting, the stochastic noise has little impact on the total number 
of inflationary efoldings even if the inflaton starts with a Planckian energy density. A “rolling” 
volume factor, i.e. weighting in proportion to the volume at that time, also leads to a monotonically 
decreasing Hubble constant and hence no eternal inflation. We show how stronger selection effects 
including a constraint on the initial and final states and weighting with the final volume factor can 
lead to a picture similar to that usually associated with eternal inflation. 
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I. INTRODUCTION 


If inflation is to provide a satisfactory explanation of the early universe, it needs both 
to hnd a successful microphysical implementation and to answer the question of why the 
universe started out in a high energy inflating state. The phenomenon of “eternal infla¬ 
tion” [1, 2] has been proposed as a solution of this second problem [3-7]. Even if the 
inflaton potential monotonically increases away from its minimum, quantum fluctuations al¬ 
low the inflaton held to jump uphill in some regions which would then expand exponentially. 
It is argued that this process, once started, can allow inhation to continue indehnitely, and 
that in all likelihood there was a great deal of inhation in the past of any observers “like 
us”. The phrase warns that selection ehects are involved, and that these might be important 
in evaluating the predictions for various inhationary models [8]. At the simplest level an 
observer “like us” could not live in a phase still undergoing high-energy inhation. So, even if 
such a phase dominates spacetime, no observers “like us” are there to see it. Attempts have 
been made to implement selection ehects with the aid of a volume-weighting, assigning more 
weight to larger regions of the universe, which are assumed to contain a greater number of 
“typical” observers. 

The anthropic leanings of the discussion can be disguised by attempting to rephrase 
selection ehects in a more physical way, for example by demanding enough inhation to 
provide superhorizon correlations on the last scattering surface say. This amounts to the 
requirement that there be something like hfty or more inhationary efolds in typical models 
(see e.g. [9]). But the observed superhorizon correlations in the universe (including large 
scale homogeneity and isotropy) appear to be much stronger than those required for suc¬ 
cessful galaxy formation. If large amounts of inhation are exponentially less likely than 
small amounts, one might interpret the existence of superhorizon correlations in our cosmic 
microwave background as evidence against inhation. 

In any case, it seems important to try and develop calculations of conditional probabilities 
within inhationary models, taking into account the back-reaction of quantum huctuations on 
the process of inhation itself. At hrst glance, one might think of a slow rolling inhaton held 
as being similar to an overdamped harmonic oscillator in the presence of weak stationary 
noise. After a short time such an oscillator ends up at the bottom of its potential and it 
only rarely huctuates appreciably upwards. Memory of initial conditions is exponentially 
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suppressed, and small fluctuations away from the minimum are exponentially more common 
than larger ones. If such a model were correct, one might expect “our” universe to only have 
the minimum possible number of efoldings consistent with the existence of a galaxy say. 

As we relate here, taking into account the held dependence of the Hubble damping and 
noise leads to a qualitatively different picture. As we shall discuss, with a change of variable 
we see that the system is actually an upside-down, over-damped harmonic oscillator, for 
which there is no stationary state. 

Nevertheless, the system can be studied via a simple Langevin equation which, for a 
particular choice of inflaton potential V = and in the slow-roll approximation, is linear 
and hence exactly soluble. Even though the dynamical evolution of the system is trivial, 
the quantities we are interested in, such as the number of inflationary efolds, are highly 
nonlinear and nonlocal in time and hence nontrivial to compute. Nevertheless, we are able 
to make progress and obtain a number of new results, extending earlier works on stochastic 
inflation [10-16]. The main advantage of the Langevin approach we take over that involving 
the Fokker-Planck equation (see e.g. [15, 17, 18]) is that quantities which are non-local in 
time, such as the number of efolds, are more easily treated analytically. We are also able to 
include various proposed weightings such as volume factors rather straighforwardly. 

Before we outline these calculations, however, let us explain the physical setup which we 
believe is approximately described by the simple stochastic model we employ. 


II. CAUSAL INTERPRETATION OF A HUBBLE VOLUME DURING INFLA¬ 
TION 

We shall be following a region whose size is the Hubble radius during inflation, and 
computing the evolution of the spatially-averaged held 0(t) in this region. The justification 
for focusing on just one such Hubble volume is that it spans the past light cone of an observer 
located far to the future. The spacetime volume inside this past light cone can be considered 
as an isolated physical system: given initial (or final, or mixed initial and final) conditions 
and a set of dynamical laws, its state should be completely describable without reference to 
the exterior. As long as causality holds it cannot be influenced by anything outside it. 

Consider the past light cone emanating from some point at time ti in a universe described 
by a flat FRW metric with scale factor a. At an earlier time t this light cone encloses a 
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FIG. 1: Heuristic conformal diagram showing an inflating patch of critical size, indicated by the 


spacetime within the cone. The physical radius of the patch at time t is the physical 

volume of the patch does not increase exponentially in time. 


sphere of physical radius 



( 1 ) 


If a{t) is increasing quasi-exponentially, i.e. the Hubble parameter H(t) = a/a is positive 
and only slowly varying with time, then the integral is dominated by its lower limit and 
rphys(t) becomes approximately equal to the Hubble radius at time t. 

So we are interested in describing the evolution of the scalar held when averaged on the 
scale of the Hubble radius at that time. We shall employ a simple stochastic model to 
describe the scalar held huctuations which is standard in discussions of eternal inhation. 
According to the model, the scalar held acquires a huctuation 6(f) ^ H on a. scale of order 
the Hubble radius, every Hubble time. We may understand the equations as describing the 
history of the huctuating held as seen in the past light cone of a point in the far future. 

Usually one hears that the physical size of a region that is known to be inhating increases 
rapidly. This statement is only accurate in the case that the initial region is so large that 
causal inhuences coming from outside the initial region cannot propagate far enough into it 
so as to shut inhation down. The critical size turns out to be the inhationary Hubble radius. 
The stochastic model applies to a region that lies on this knife edge: small enough to he 
within the past light cone of a future observer, and large enough to remain in an inhating 
state well into the future. See Fig. 1. 
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III. BACKGROUND EVOLUTION 


Our starting point for the stochastic approach to inflation will be the Friedmann and 
scalar held equations for a hat FRW metric: 




8nG 


1 .2 


<p + 


0 , 

0 . 


( 2 ) 

( 3 ) 


One then considers linearized perturbations about a background solution. In each Hubble 
time, new quantum huctuations in the scalar held are generated and freeze out on the scale 
of the Hubble radius with an amplitude of order H [2], 

We model the ehects of such huctuations by adding a stochastic noise term onto the 
right-hand-side (RHS) of Eq. (3). We take the noise to be proportional to delta-function- 
normalized Gaussian white noise, n{t), which obeys 

{n{t)) = 0, {n{t)n{t')) = 5{t -1'). (4) 

Clearly, n{t) has dimensions of mass to the power one half. The coefficient may be deter¬ 
mined, up to a numerical coefficient of order unity, by dimensional analysis. The only scale 
entering into the huctuations is H, with dimensions of mass. Thus we need on the 

RHS of (3) for the correct dimensions. We shall not be concerned about the numerical coef- 
hcient here, but if desired this can be determined by normalizing the model to a calculated 
quantum correlation function (see e.g. [10]). 

There is a long history of attempts to to take account of the ehects of the huctuations 
on the background evolution in inhation, going back to [19] for example. Ref. [2] develops a 
picture of the held evolution over a hat region of its potential as Brownian motion. It was 
later realized that this motion could also be important for helds with unbounded potentials 
and that this ehect might be crucial in predicting what an observer might expect to observe 
(see [3-5]). Numerical simulations of this process have been performed (see e.g. [16-18]). 

One can often self-consistently make the slow-roll approximation, even in the presence of 
noise, which involves neglecting the 0^/2 term in (2) and the 0 term in the noisy version 
of (3), leaving us with 

3H<j) + V^ = (5) 
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where H = H{(j)) = A/87rGV^(0)/3. From now on we shall adopt reduced Planck units, 
setting 8nG = 1. 

IV. EXPLICIT SOLUTION FOR SLOW-ROLL \cf)^ 

There is a special choice of potential for which (5) simplihes, namely for V = as 
noticed by Nambu [13]. This is an interesting potential to investigate in its own right and 
might reasonably be expected to be representative of other models with simple power law 
potentials such as After a change of variable, the Langevin equation becomes linear, 

with held-independent noise. Remarkably, the new stochastic variable is just the physical 
Hubble radius. 

Dropping the second-derivative term, dehning R = l/(\/A/30^), the physical Hubble 
radius in the slow-roll approximation, and introducing a = 8a/A/ 3 and (3 = 2^A/3/3 we 
have 

R — aR = (3n{t), (6) 

where we have changed the sign of the RHS relative to (5). As introduced above, n{t) is 
delta-normalized Gaussian white noise with zero mean. We use angular brackets ((...)) to 
denote the ensemble average of a quantity (...) over histories of the noise function. We 
shall use double angle brackets (((...)'^)) to denote the gth cumulant of the distribution of a 
quantity (...), and also for connected correlation functions of products of quantities. Note 
that powers of f3 count the number of times that the noise enters into any expression. 

Eq. (6) describes, as claimed, an over-damped, upside-down harmonic oscillator with 
linear noise. It immediately reveals a potential problem in that nothing prevents R from 
crossing zero for some trajectories. Both the stochastic model and the underlying held theory 
of inhation may be expected to break down there, since the scalar held tends to inhnity. A 
similar pathology occurs in the Fokker-Planck approach [13, 15], and we can try to apply 
analogous workarounds here where necessary. 

If we specify an initial condition i? = rg at t = 0 say, we have the integral solution 

R{t) = roe"* (^1 + ^ J . (7) 
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Averaging over the noise, we find 


{R{t)) = roe“‘, 
{R{t)R{t')) = 





2a min(t, 



( 8 ) 

( 9 ) 


For the second expectation value we have multiplied two integral solutions together, taken 
the ensemble average using the assumed properties of the noise, given in (4), and then 
performed the integrals over the dummy time variables. 

The mean value {R{t)) does not involve the noise at all. Changing variables back to (j) 
we see this just represents the classical slow-roll solution 0 = . Because the noise 

is Gaussian, the /^-distribution is also Gaussian, with a mean p = roe“* from (8) and a 
variance — 1)/(2 q;) given by setting t' = tin (9). More generally, one can derive 

a simple integral expression for the expectation value of any function f{R{t)) of R that can 
be expressed as a Fourier integral, by taking the averaging inside the Fourier integral and 
writing in terms of cumulants as _ 

Defining as a dimensionless measure, of the variance in R{t)\ 


X\t) ^ 


{{R\t))) 

{R{t)f 


we see that R = 0 lies 1/X standard deviations to the left of {R{t)). Thus if X is small, 
one can hope that for reasonable weightings the problems associated with R = 0 can be 
neglected. X is small at small times for all starting values ro, and remains small for all time 
if ro is large enough that X^ = X^(cx)) = P‘^/{2arl) 1. is a quantity with a nice 

physical meaning: it is, up to a numerical coefficient, the ratio of the initial energy density 
to the Planck energy density. So if we start well below the Planck scale, we do not expect 
subtleties at i? = 0 to be important when following a patch forward in time. 

Strictly speaking, one should not trust the theory at all if the density approaches the 
Planck density. Nevertheless, one can attempt to patch it up by imposing reflecting [13] or 
absorbing [15] boundary conditions at i? = 0. We shall discuss the latter in Sec. VI below. 
In our Langevin treatment these two prescriptions correspond respectively to averaging or 
differencing the probability densities over trajectories that start at both tq and —rg. 
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V. ASYMPTOTIC EXPANSION FOR THE NUMBER OF EFOLDINGS 


Let us use our integral solution to calculate the expected number of inflationary efolds 
N{T) = dtH{t). If we were following a comoving region of a homogeneous inflating 
universe rather than one physical Hubble volume in an inhomogeneous universe, N{T) would 
simply be the number of inflationary efolds that the region had experienced up to time T. 
We shall continue to call N the number of efolds, keeping in mind though that it bears 
no simple relation to the physical size of the volume that we are actually tracking in the 
calculation. 

Since H is just 1/i? in our approximation, 

(N[T))=j'^dt =J^dt (1 + /(«))■') . (11) 

where 


/(i) ^ ^ f 

ro Jo 


dtie 


( 12 ) 


We now perform a formal Taylor expansion in I{t) of the reciprocal inside (11). The justi- 
hcation for this procedure is that as long as the noise term can be treated as “small”, an 
expansion in the noise makes sense. 

We can take the expectation value of powers of I term by term making use of the fact 
that the noise is assumed to be Gaussian so that all higher-order correlation functions of it 
can be expressed in terms of products of its two-point function via Wick’s theorem. It thus 
turns out that (/'^(t)) vanishes if q is odd, and equals (2p— l)!!X^^(t) if g = 2p is even, with 
X{t) as dehned above and {2p — 1)!! = {2p — l){2p — 3)... 1. Thus 


^—at 


{H(,t)) = — (1 + X^t) + SX-'O) + ... + (2p - 1)!!X2!'(() + ...). 


(13) 


Finally, we can perform the integral over t term-by-term to obtain a series for (iV(T)). 

This procedure works particularly neatly if we take t oo, since then the time integrals 
are straightforward, using dye~^(l — = 2Pp\/{2p+ 1)!!. With = -/V(cx)), one 

obtains: 


(lVo„) = lVd(l + l'U’fi + 3.^.Ay + ... + ^^Ayf+,,,^, (14) 

where Wi = ^ is the “classical” number of efolds expected in the slow-roll approximation 
at late times in the absence of noise. 



Unfortunately the two series expansions (13) and (14) diverge. We shall investigate this 
divergence in more detail in the next section. For now, let us understand these series as 
asymptotic expansions in terms of the initial energy density. This is plausible since the 
Taylor expansion of (1 + inside the integral cannot be valid for all histories of the 

noise. But if the noise is weak enough we do not expect this inaccuracy to be signihcant. 
An analogous example is the calculation of the expectation value of (1 + ex‘^)~^ or (1 — ex)~^ 
say for a Gaussian distribution with zero mean and unit variance in the case that e -C 1, for 
which the perturbative expansion, although only asymptotic, is still extremely accurate. 

Our second illustration (1 — ex)~^ above is instructive in suggesting a formal procedure 
to account for effects associated with i? = 0. Here we have a divergence in the calculation of 
the expectation value, at a; = —1/e, far into the tail of the distribution. If we don’t trust the 
detailed form of the distribution for large negative x, we may feel that such a divergence is 
unphysical and irrelevant. A natural mathematical way of removing such a divergence, and 
allowing us nevertheless to continue to use the Gaussian form, is to calculate the expectation 
value of principal values (denoted P.V.) of linearly-divergent quantities rather than the 
quantities themselves. In the stochastic model this prescription turns out to be related to 
imposing absorbing boundary conditions at i? = 0. 

After these comments, let us see what we can glean from Eq. (14). We see that fluctu¬ 
ations increase the expected number of efolds that a patch undergoes, but only by a small 
multiplicative fraction if the energy density starts low compared to the Planck density. Sig¬ 
nihcant corrections to the classical result only occur at such high energy densities that we 
can have no conhdence in the applicability of the theory. 


VI. THE TOTAL NUMBER OF EFOLDINGS 

We would like to investigate in greater depth the divergence of our series for {H(t)) and 
(iV(T)). In the Fokker-Planck approach, a natural way to regulate the problem at i? = 0 is 
to impose the “absorbing” boundary condition, that the probability is zero at i? = 0 [15]. 
This can be done for a Gaussian probability distribution by rehecting it around R = 0 and 
subtracting the rehected copy off the original one. This new probability distribution must 
then be renormalized. In the Langevin approach this procedure corresponds to ignoring any 
paths that reach R = 0. One could alternatively impose “reflecting” boundary conditions for 
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which the gradient of the probability distribution at i? = 0 is zero, or indeed consider a mix¬ 
ture of the two conditions. Only the pure absorbing condition ensures that the probability 
vanishes at i? = 0, and for definiteness we shall focus on this case here. 

From our Langevin approach, we know that before the boundary condition is imposed 
R is Gaussian-distributed, with mean fi and variance given above. Applying the reflect- 
and-renormalize procedure outlined above, the expectation value of some function f{R) is 
given by 


J^dxfix) (e 


(/(^))abs = 


/o+ dx 


e 2 <t^ — e 


(The subscript “abs” is for absorbing.) We extend the range of / from the positive reals to 
the entire real line by demanding that / is odd. Then we let x —x in the second term 
of the upper integral, and divide top and bottom by The numerator is then the 

expectation value of the P.V. of / with respect to the Gaussian we started with over the 
entire real line. The denominator is the error function eri{jj,/{\/2a)) = erf(l/(v^W(t))). 

Since H(t) = l/R(t) we can apply this to work out (^f(^))abs- 



To make contact with our asymptotic expansion above, we write P.V. s-s a Fourier 

integral and then take its Gaussian expectation value as discussed earlier. This gives the 
numerator as the integral 

(p-''- )=^ ■ VM r "" (^) 

If we then continually integrate by parts on the sin we generate the asymptotic expan¬ 
sion (13). 

We still have to consider the denominator. Now erf(l/(v^X(f)) = 1 — erfc(l/(v^X(t))), 
and for small X(t) erfc(l/(v^X(t)) ~ a/2/ ttX (t) d)) (see e.g [20]) which is exponen¬ 
tially small. So the denominator itself is exponentially close to one, and the entire expression 
is accurately approximated by the asymptotic expansion (13). 

Thus we see that effects coming from imposing boundary conditions at i? = 0 are negli¬ 
gible to the extent that X{t) is small. 
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This condition at i? = 0 does allow us to go on to investigate the large X regime also. 
We may rewrite (15), applied to ff(t) = 1/R(t), as 


where 


(mu. = 


r f 

/“ dxW(x) 


(18) 


W{x) = sinh e 


(19) 


By rescaling the integration variable, one hnds that the expected Hubble constant can be 
expressed in the form 

(H(t)).,.. = (20) 

One hnds the asymptotic behaviour: 

F(z) ~ 2 :, for z ^ 1; (21) 

F(z) ~ \/(7rz)l2, for z :$> 1. (22) 


Using this, one estimates the expected total number of efoldings (IVoo)abs ~ /o°° (-^(^))abs 

as: 

(^oo)abs ~ ^ 

CXi Q O 

which is just the standard classical result, and 

/ (52 ^-1/2 / 7 r\ 3/2 1 

( 2 ) forr„«l, (24) 

which up to a numerical constant is just the standard classical slow-roll result for 

inhation starting at the Planck density. Hence we conclude that quantum huctuations, 
treated as stochastic noise, do not qualitatively alter the expected number of inhationary 
efoldings beyond the classical result, for initial densities right up to the Planck density. 


VII. “ETERNAL INFLATION” FROM STRONG SELECTION 

Imagine we know that at two times the scalar held was well up the hill within the Hubble 
volume we are tracking. What can we say about the likely value of the held at intermediate 
times? Did the scalar held roll downhill and then huctuate back up, or did it huctuate up 
and then roll down to its prescribed hnal value? 
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There are two equivalent ways of calculating such an effect. One is to use Bayes’ theorem 
explicitly. The other to exploit the fact that the noise is Gaussian, which means that the 
full probability distribution for multiple events is given by exponentiating their covariance 
matrix (this method cannot be used if the boundary conditions at i? = 0 are important). 
We shall present the former method here. 

Let us take 0 < t < T. We denote by pi|i(rt|ro)(ir the probability that at time t, R lies 
between rt and rt + dr given that i?(0) = ro; and by Pi\ 2 (rt\rT, ro)dr the probability that at 
time f, R lies between and + dr given that R{T) = r^p and i?(0) = tq. 

Then, by Bayes’ theorem. 


Pi|2(n|rr,ro) = 


Pi\iirt\ro)Pi\2irT\rt, ro) 


(25) 


Pi\i{rT\ro) 

and we know what all the terms on the RHS are, since Pi\ 2 {rT\rt,ro) = Pi|i(rr|rt). 

If we neglect R = 0 effects, the probability distributions are Gaussian and we can use our 
formulae (8) and (9) to evaluate the means and variances for the RHS. With the constraints 
on i? at 0 and T we hnd that at the intermediate time t, R is again Gaussian-distributed 
with a t-dependent mean of 


{«(«))„„ = 


ro sinh q;(T — t) + sinh at 


(26) 


sinh aT 

(The subscript “con” is for constrained.) In the Appendix we discuss an interesting formal 
property of {R{t))^^^, which is that it obeys a second order differential equation related to 
the original hrst order equation (6) in a simple way. 

One can also easily hnd the variance in R by looking at the coefficient of the term 
quadratic in in the exponent on the RHS of (25). It turns out to be: 

2 sinh a{T — t) sinh at 


mm... = 


(27) 


' 2q; sinh aT 

Interestingly, this is independent of r-Q and r-p, and attains its maximum value of 
(/5^/2 q;) tanh(Q;T/2) at the midpoint t = T/2. This is less than P'^/2a, so the standard 
deviation in R here is always bonnded by a constant of order nnity. As long as {R(t))^^^ is 
always signihcantly larger than unity, we expect that we may neglect effects from R = 0. 

Furthermore, if the minimum of {R{t))^^^ is much less than r^, then relatively close to T 
{R{t))^.^^ behaves like independent of tq. This is the same as the noise-free slow- 

roll solution which passes throngh r-p at time T, and which in the far past would approach 
arbitrarily small R and thus large 0. 
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FIG. 2; A plot showing the expectation of ii = 1 / {y\/^(jP‘) against time for paths constrained to 
take specified values at two different times. In the example shown, R = 8 at at = 0 and at = 6 
(solid line). Also indicated are one and two sigma contours for the trajectories (dotted lines), and 
the noise-free solution for R that passes through R = 8 at at = Q (dashed line). We have set 
0^/2a = 1/36 in making this plot. 

For example, if we take = ro, {Rif)) con takes a minimum value of ro/ cosh(Q;T/2) at 
t = T/2. Demanding that the minimum value is greater than the standard deviation at that 
time puts a condition on the values of T for which we may reliably neglect effects associated 
with the i? = 0 boundary: this condition reads sinhaT < 2/X^. Since is typically 
small if the held starts low down the potential, we see that as long as T < ln(4/X^)/Q; 
boundary effects from i? = 0 should be unimportant. As long as this is so, we can now 
answer the question posed at the beginning of this section. If the held is known to be at 
the same place on the potential at times 0 and T, the Hubble radius lies at smaller values 
in the interim, corresponding to the held being further up the potential. See Fig. 2. 

So we have shown that if inhation has lasted a reasonably long time (i.e. we know that 
the held was displaced well up its potential at two widely separated times), the held was 
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likely to have rolled up the hill to higher held values, and then turned round and rolled 
back down, approaching the standard slow-roll solution. This behaviour shares many of the 
characteristics that are commonly associated with the phrase “eternal inhation”. 

As an example, imagine we select paths in which the scalar held takes a value of order 
the Planck mass at two widely separated times. Recall that this is the value for which the 
slow-roll approximation fails, near the end of inhation. The above calculation shows that 
the expected trajectory is one in which the held hrst runs uphill to very large values before 
rolling down along a nearly classical slow-roll path to the hnal value. For such paths, one 
would expect the usual predictions for the observable density perturbations. 

However, by demanding that T is large we have imposed a very strong selection of paths; 
the vast majority of paths run downhill. Now let us consider weakening the selection by 
reducing the time T between our initial and hnal states. Even if we insist upon hfty efoldings 
of inhation, the likely trajectory is very diherent from the classical slow-roll one usually 
considered. In standard inhation where one ignores the ehect of quantum huctuations on 
the background, the density perturbation amplitude is approximately However, in 

the selected ensemble the velocity 0 of the scalar held vanishes at an intermediate time. This 
would lead one to suspect that for the selected ensemble, density perturbations which exited 
the Hubble radius around the time T/2 would be far larger than those in standard inhation. 
Unless T is chosen to give substantially more than hfty efoldings, the normal predictions of 
inhationary perturbations will not be obtained. 

To conclude this section, we have seen how selection ehects can indeed make a scalar 
held roll uphill before rolling roll back along a near-classical slow-roll path. This is the same 
type of behaviour as that sought from eternal inhation, but we have obtained it through 
an imposed ad-hoc selection, rather than from volume weighting, the ehect to which it is 
normally attributed. We now turn to an investigation of volume weighting to see whether 
this can produce a similar ehect. 

VIII. ETERNAL INFLATION AS A VOLUME EFFECT 

The usual argument for inhation being eternal to the future with an inhaton held with an 
unbounded potential is summarized and criticized in [21]. The main idea is that regions of 
space in which the held huctuates higher inhate more rapidly, and thus after a long amount 
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of time the majority of the physical volume is dominated by inflating regions at the Planck 
density. Occasional regions fluctuate downwards and inflation stops, and we are expected 
to live in one of these regions. The main criticism is that this line of argument is not gauge 
invariant, i.e. comparisons of volumes depend on the choice of time-slicing. 

We argued above that the stochastic equation really only follows the evolution of the field 
inside one physical Hubble volume. The common supposition is that it rather describes the 
average behaviour of the held in a hxed comoving volume. In Fig. 1, this corresponds to 
assuming that the held on the constant-time slice through the entire conformal cylinder is 
the same as it is on the piece of the slice inside the cone. In this case, one may then argue 
that expectation values of various quantities evaluated at time t should be calculated by 
weighting each trajectory by its volume factor In this section, we shall explore the 

consequences of this suggestion. 

We thus dehne the “volume-weighted” average of some quantity (...) at time t via: 

• ( 28 ) 

(The reasoning behind the subscript will become clear shortly.) 

Let us attempt to calculate the volume-weighted number of efolds This is an 

interesting example because the usual claim (see e.g. [6]) is that if 0o is above a critical value 
of order then the physical volume of space that is inflating should increase withont 

limit, corresponding to eternal inflation. Thus we would expect to see a qualitative change 
in the behaviour of around this critical starting valne. 

We start by introducing the formal generating fnnction where jj, 

is a nnmber. Differentiating this with respect to fi and then setting p = 3 gives us the 
volume-weighted number of efolds, (hence the snbscript in (28)). We may think 

of IFt(/i) as 

CXD 12 

^^((iV(i)'.)), (29) 

1 

n=l 

So we can express (iV(t))^^ 3 j in terms of the connected moments or cnmulants of N. The 
nth cumulant of N can be obtained as usual from the regular moments of N, which can be 
calculated order by order in the noise. (Such a procedure is described in more detail in the 
following section where it is needed for a wider-ranging calculation.) 
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In the limit t oo, it turns out that the nth cumulant goes like times a power series 
in starting at the (n — l)th power of X^. So to order X^ we hnd 

{{N^)) = /Vd (l + \xl + ^ j (30) 

({Nl)) = Nl + 3Ai + . ^ (31) 

{{Nl)) = Nl + ,,,) (32) 

with all higher cumulants higher order in X"^. Following the aforementioned procedure and 
collecting terms in we hnd 



Thus we obtain the fractional correction to the classical number of efolds Nc\ to fourth order 
in the noise. One sees that for Nc\ 3> 1 (but still 1), we in fact seem to have an 

expansion in X^X^ rather than in just X^. So our formula gives only a small fractional 
correction to the number of efolds as long as XdX^ 1. 

But XciX^ = {2a‘^rQ) ~ A0o and so we see that when 00 ~ the volume-weighted 

number of efolds begins to change substantially from the noise-free result. This is just the 
criterion mentioned above that is usually given for the onset of eternal inhation. 

So our method has established some contact with the usual approach to eternal inflation. 
However, we caution that the assumption that one should volume-weight is not yet (in our 
view) well-founded, since the dynamical equations being used are only following the value 
of the field in one Hubble volume. 

IX. A CLOSER INVESTIGATION OF VOLUME-WEIGHTING 

In the above section we introduced volume-weighted averages of quantities dehned at 
a given time, where the volume-weighting was performed at that time. We used this to 
calculate the volume-weighted number of efolds at a given time f, and then considered 
taking the limit t ^ oo. 

In this section, we shall consider two sorts of volume weighting. If one believes in volume¬ 
weighting as a physical spatial averaging, then to compute observables at time t one might 
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want to weight with the volume at that time. On the other hand, if one thinks of volume¬ 
weighting as a selection effect for “typical” observers, one might want to weight with the 
“hnal” volume, at some much later time. To cover both cases, we shall compute the Hubble 
parameter as a function of time, H(ta), weighted by the volume at a later time t > ta- 
Taking the limit ta —^ t, we obtain the trajectory of the Hubble parameter with a rolling 
volume weighting. If instead t is taken to cx), we obtain the expected history of the Hubble 
parameter when selected by hnal volume. This gives what a “typical” late-time volume- 
weighted observer should expect to see in his/her past. (Note that our treatment in the 
previous section took the “spatial averaging” point of view but focused only on the hnal 
number of efoldings.) 

To calculate such averages involving the Hubble parameter, one needs to generalize the 
generating function introduced above, promoting p from a number to a function of time. 
Thus one dehnes the generating functional 

H^t[p] = ln(e'^‘[^]>, (34) 

where 

Mtln] = f dtifi{ti)H{ti). (35) 

Jo 

Now by functionally diherentiating with respect to and then setting p = 3 for all time, 

one obtains volume-weighted cumulants involving H{ta)- Diherent choices for p correspond 
to more general forms of weighting, and we may denote expectation values with respect to 
them by ((.. t should we wish to use them. Note that p = 0 corresponds to the natural 

weighting. It can be useful to take p to be an arbitrary constant at intermediate stages of 
the calculation, and only set p = 3 at the very end. Then one can easily identify any ehects 
coming from the volume weighting. 

In order to proceed perturbatively in the noise, we again expand W in terms of the regular 
(i.e. non volume-weighted) cumulants of M[p]: 

OO 

= E ;;t ■ (36) 

n=l 

These cumulants involve multiple integrals J dti... J dtj over products of the Hs. Each H 
is written as: 

Hit.) =- (1 -/, + /? - If + (37) 

ro 
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where li is a shorthand for I{ti) as defined in (12). We use Wick’s theorem to express 
higher-order moments of I in terms of the two-point function (/ 1 / 2 ), which is given by 

(/ 1 / 2 ) = (1 - _ (38) 

Each power of Xqo counts an order of the noise. As in the previous section, it turns out that 
the nth cumulant is of order 2(n — 1) in the noise. 

Let us now use the above machinery to calculate {H(ta))^ to second order in the noise. 
We need the first and second cumulants of M, each evaluated to second order in the noise. 
Functionally differentiating the former leaves us with {H{ta)), or 

-(l + X^(l-e-2“‘“)). (39) 

^0 

Functionally differentiating the latter cumulant yields 2 dtin{ti) {{H{ti)H{ta))), or 

rt p-cdl p-ata 

2 / dtifi{ti) - (1 - _ (4Q) 

Jo ro ro ^ 

This can be evaluated assuming p to be constant in time. Utilizing these results and setting 
/i = 3 leads to 

p-ata ( 

+ 3XdX^ (^1 - e-“*-^-(e""*" - e”"*) j \. (41) 

As discussed above, we are interested in two cases. The hrst, involves taking ta —> t. This 
choice implements a “rolling volume weighting”, where at any time the volume-weighted 
average is performed using the volume at that time. This choice would correspond to what 
one might usually consider to be a spatial average over physical volume. We obtain 

({ff(i))>,..3,, = + -’ft (1 - e-“) + 3NaXl (^1 - e"" - |. (42) 

(This result can also be obtained by calculating (X(t))^^ 3 j as discussed in the previous 
section to second order in the noise and then differentiating with respect to t.) A Taylor 
expansion of the correction term coming from the volume-weighting starts at 0{{atY). This 
is to be compared with the 0{at) variation coming from the expansion of the “classical 
rolling” term. One (initial) Hubble time is tn = ^/Hq = r^, so atu = OiVQ = l/Xd, which is 
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small if the classical number of efoldings is large. Thus as long as one is interested in times up 
to of order a few Hubble times from the start, rolling volume-weighting does very little even 
if iVciX^ ~ A0 q 3> 1. One only starts getting corrections of order when at ~ 1, that 

is at times for which, without noise, a substantial fraction of the total number Nc\ of efolds 
would have occurred. We have checked that including the fourth order term does not affect 
this conclusion, and strongly suspect that neither will further higher order terms. Hence 
we conclude that a rolling volume weighting leads to a monotonically decreasing Hubble 
constant, and will not produce the uphill motion required by the eternal inflation picture. 

Now let us consider the second limit. This is to take t oo. Thus we are weighting 
trajectories by their dual volume. As discussed above, this in some sense corresponds to 
the history of the held that a physical observer at late times might expect to see in his/her 
past light cone, assuming that physical observers are evenly distributed over the final three- 
volume. We obtain: 

= —{l + (1 - e-““) + 2/VdA'^ (1 - e-“-) }, (43) 

Now the volume-weighting term has a Taylor expansion that starts at 0{ata), the same 
order as that from the “classical” exponential. So now for large enough Wi, we can have 
the expectation value of the Hubble parameter increase with time. The requirement is that 
QNciX^ > 1, or 00 ^ again the criterion for eternal inflation. By ata ~ 1, the 

fractional correction to the noise-free result is becoming large and we do not trust our result 
in the details. Of course, if the noise is very important we should not overly trust the t ^ oo 
limit in the first place. 

We can investigate the spread in the final-volume-weighted trajectories by considering 
the ratio The lowest-order term in the numerator is already 

O(X^), so we must calculate the numerator to fourth order to obtain the fractional change 
in the ratio due to volume-averaging to second order. Thus we need the third cumulant of 
H for the numerator, but the calculation can be performed. Taking ata —^ cxo to obtain the 


late-time behaviour and dropping terms down by Wi we find 

(1 + y Ady (44) 

{A(oo)>„.3,oo V 5 ) 

Thus we see that in the regime of eternal inflation the noise is playing an important role. 
One might be tempted to conclude that for 0o ~ A“^/® the volume-weighted system becomes 


“noise-dominated” and a “non-perturbative” way of treating the noise is called for. 
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X. CONCLUSIONS 


We have seen that inflation with quantum jumps may be modelled as an upside-down, 
overdamped harmonic oscillator with noise. The critical role that selection effects have on 
expectations for the field’s observed history has been clearly demonstrated. In the absence 
of any selection effects, we have seen that the presence of stochastic noise hardly alters the 
motion of the inflaton held downhill; the total number of inflationary efolds is essentially 
the same as the classical result. Through strong selection one can obtain trajectories which 
travel uphill; if the held is constrained to he at a particular inflating value at two widely 
separated times we have shown that it is likely in the interim to be at higher field values. 

We also investigated volume-weighting in our approach. With a rolling weighting, such 
as might be attributed to averaging over physical volume, we And no significant effect on 
the Hubble parameter’s evolution until late times, for which inflation would classically have 
completed a substantial fraction of its total number of efolds. If, alternatively, histories 
are weighted by the final volume, qualitative changes in the mean trajectory of the system 
can occur within a few Hubble times of the start. However in this situation, it is not clear 
whether the perturbative treatment is valid. 

Our techniques might also be applicable for more general forms of averaging. For ex¬ 
ample, one might consider the possibility of weighting by 4-volume. This would be the 
right thing to do if observers spring into being from inflating regions in a stochastic manner 
consistent with Poisson statistics (say by bubble nucleation), and that the appearance of 
a region supporting physical observers has a negligible effect on the background evolution 
(unlike bubble nucleation). If the 4-volume were finite, the result would be gauge-invariant 
and ambiguity-free also. The generating function could then be In and dif¬ 

ferentiating this with respect to ju and then setting /i = 3 would give the 4-volume-weighted 
number of efolds. The generating function could be calculated order-by-order in the noise 
using our methods. As we have explained, however, we are not convinced that any particular 
form of volume-weighting is physically correct, since the stochastic equation being used only 
describes the time evolution of a single physical Hubble volume. 

There are several possible extensions of the work reported here. Numerical simulations 
could be used to extend the perturbative treatment we have given and to check the generality 
of the behaviour found here for other potentials. It would be interesting to know whether 
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inflation models with much flatter potentials exhibit the same qualitative behaviour. It 
would also be interesting to extend the model to allow the inflaton held to jump back up 
from its minimum into an inhating regime, a situation referred to by Garriga and Vilenkin as 
a “recycling universe” [22] (see also [23]). In this situation one might hope to hnd a steady 
state solution, although the typical bout of inhation would involve very small numbers of 
efolds and so such a model would be observationally disfavoured. 
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APPENDIX; A PROPERTY OF 

In this Appendix we discuss an interesting fact about the constrained expectation 
which is that it satishes the second order deterministic differential equation 

where the operator appearing is the product of the original operator in the equation of 
motion, (6), and its time reverse. That this should be so is in fact a general consequence 
of time-translation invariance of the operator appearing in the original equation of motion, 
(6), and of the noise ensemble. The point is that the constrained expectation is 

proportional to the correlator {R{t)R{T)), as may be seen by writing the joint (Gaussian) 
probability distribution for R(t) and R(T). To show the correlator obeys the stated equation 
(A.l), note first that 

for all t ^ T, since from the equation of motion (6) this equals {n{t)n{T)), which is 
zero at unequal times. Now, we can write the left hand side as (^ — a) acting upon 
P{n{t)R(T)). But R(T) is a sum of a particular integral which does not correlate with 
the noise, and an integral (3 dt'G{T,t')n{t'), where G{T,t') is the Green function of the 

original operator. So {n(t)R(T)) is just PG(T,t). Due to time translation invariance, the 
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latter is a function only of the time difference t — T. Therefore, we may replace ^ in (A.2) 
with — obtaining (A.l). The arguments we have used are rather general: it is clear for 
example that constrained expectation values related to a second order stochastic differential 
equation would obey a fourth order deterministic differential equation. 
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